Non-monotonic, distance-dependent relaxation of water in reverse micelles: Propagation of surface induced frustration along hydrogen bond networks

Layer-wise, distance-dependent orientational relaxation of water confined in reverse micelles (RM) is studied using theoretical and computational tools. We use both a newly constructed ``spins on a ring'' (SOR) Ising-type model (with Shore-Zwanzig rotational dynamics) and atomistic simulations with explicit water. Our study explores the effect of reverse micelle size and role of intermolecular correlations, compromised by the presence of a highly polar surface, on the distance (from the interface) dependence of water relaxation. The ``spins on a ring'' model can capture some aspects of distance dependence of relaxation, such as acceleration of orientational relaxation at intermediate layers. In atomistic simulations, layer-wise decomposition of hydrogen bond formation pattern clearly reveals that hydrogen bond arrangement of water at a certain distance away from the surface can remain frustrated due to the interaction with the polar surface head groups. This layer-wise analysis also reveals the presence of a non-monotonic slow relaxation component which can be attributed to this frustration effect and which is accentuated in small to intermediate size RMs. For large size RMs, the long time component decreases monotonically from the interface to the interior of the RMs with slowest relaxation observed at the interface. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4732095]

Physical characterization of the rack effect and hydrogen bond networks in blue copper proteins

A summary of previous research is presented that indicates that the purpose of a blue copper protein's fold and hydrogen bond network, aka, the rack effect, enforce a copper(II) geometry around the copper(I) ion in the metal site. In several blue copper proteins, the C-terminal histidine ligand becomes protonated and detaches from the copper in the reduced forms. Mutants of amicyanin from Paracoccus denitrificans were made to alter the hydrogen bond network and quantify the rack effect by pK_a shifts.

The pK_a's of mutant amicyanins have been measured by pH-dependent electrochemistry. P94F and P94A mutations loosen the Northern loop, allowing the reduced copper to adopt a relaxed conformation: the ability to relax drives the reduction potentials up. The measured potentials are 265 (wild type), 380 (P94A), and 415 (P94F) mV vs. NHE. The measured pK_a's are 7.0 (wild type), 6.3 (P94A), and 5.0 (P94F). The additional hydrogen bond to the thiolate in the mutants is indicated by a red-shift in the blue copper absorption and an increase in the parallel hyperfine splitting in the EPR spectrum. This hydrogen bond is invoked as the cause for the increased stability of the C-terminal imidazole.

Melting curves give a measure of the thermal stability of the protein. A thermodynamic intermediate with pH-dependent reversibility is revealed. Comparisons with the electrochemistry and apoamicyanin suggest that the intermediate involves the region of the protein near the metal site. This region is destabilized in the P94F mutant; coupled with the evidence that the imidazole is stabilized under the same conditions confirms an original concept of the rack effect: a high energy configuration is stabilized at a cost to the rest of the protein.

Mastering the canonical loop of serine protease inhibitors : enhancing potency by optimising the internal hydrogen bond network

Background Canonical serine protease inhibitors commonly bind to their targets through a rigid loop stabilised by an internal hydrogen bond network and disulfide bond(s). The smallest of these is sunflower trypsin inhibitor (SFTI-1), a potent and broad-range protease inhibitor. Recently, we re-engineered the contact β-sheet of SFTI-1 to produce a selective inhibitor of kallikrein-related peptidase 4 (KLK4), a protease associated with prostate cancer progression. However, modifications in the binding loop to achieve specificity may compromise structural rigidity and prevent re-engineered inhibitors from reaching optimal binding affinity. Methodology/Principal Findings In this study, the effect of amino acid substitutions on the internal hydrogen bonding network of SFTI were investigated using an in silico screen of inhibitor variants in complex with KLK4 or trypsin. Substitutions favouring internal hydrogen bond formation directly correlated with increased potency of inhibition in vitro. This produced a second generation inhibitor (SFTI-FCQR Asn14) which displayed both a 125-fold increased capacity to inhibit KLK4 (Ki = 0.0386±0.0060 nM) and enhanced selectivity over off-target serine proteases. Further, SFTI-FCQR Asn14 was stable in cell culture and bioavailable in mice when administered by intraperitoneal perfusion. Conclusion/Significance These findings highlight the importance of conserving structural rigidity of the binding loop in addition to optimising protease/inhibitor contacts when re-engineering canonical serine protease inhibitors.

Vibration spectra of ferroelectric alkali trihydrogen selenites and the nature of the hydrogen bond potential in these crystals

Raman spectra of the ferroelectric LiH3 (SeO3)2 and NaH3(SeO3)2 and the anti-ferroelectric KH3 (SeO3)2 have been recorded at room temperature using a He-Ne and also an Ar-ion laser source. The infrared absorption spectra of these crystals and their deuterated analogues have been recorded in the region 400–4000 cm−1 both below and above the Curie temperature. From an analysis of the spectrum in the region 400–900 cm−1 it is concluded that (i) in LiH3 (SeO3)2 the protons are ordered in an asymmetric double minimum potential with a low barrier and the spectrum can be interpreted in terms of HSeO3− and H2SeO3 vibrations, (ii) in NaH3 (SeO3)2 all three protons occupy a single minimum potential at room temperature and below the transition temperature the groups HSeO3− and H2SeO3 are present, (iii) the proton at the inversion centre in KH3(SeO3)2 is in a broad troughed potential well and the low temperature spectrum is more likely to be due to H3SeO3+ and SeO32− species. This deviation of the spectrum from that of the previous two crystals is attributed to the difference in H-bond scheme and hence the absence of any cooperative motion of protons in this crystal.

The hydrogen bond: a molecular beam microwave spectroscopist's view with a universal appeal

In this manuscript, we propose a criterion for a weakly bound complex formed in a supersonic beam to be characterized as a `hydrogen bonded complex'. For a `hydrogen bonded complex', the zero point energy along any large amplitude vibrational coordinate that destroys the orientational preference for the hydrogen bond should be significantly below the barrier along that coordinate so that there is at least one bound level. These are vibrational modes that do not lead to the breakdown of the complex as a whole. If the zero point level is higher than the barrier, the `hydrogen bond' would not be able to stabilize the orientation which favors it and it is no longer sensible to characterize a complex as hydrogen bonded. Four complexes, Ar-2-H2O, Ar-2-H2S, C2H4-H2O and C2H4-H2S, were chosen for investigations. Zero point energies and barriers for large amplitude motions were calculated at a reasonable level of calculation, MP2(full)/aug-cc-pVTZ, for all these complexes. Atoms in molecules (AIM) theoretical analyses of these complexes were carried out as well. All these complexes would be considered hydrogen bonded according to the AIM theoretical criteria suggested by Koch and Popelier for C-H center dot center dot center dot O hydrogen bonds (U. Koch and P. L. A. Popelier, J. Phys. Chem., 1995, 99, 9747), which has been widely and, at times, incorrectly used for all types of contacts involving H. It is shown that, according to the criterion proposed here, the Ar-2-H2O/H2S complexes are not hydrogen bonded even at zero kelvin and C2H4-H2O/H2S complexes are. This analysis can naturally be extended to all temperatures. It can explain the recent experimental observations on crystal structures of H2S at various conditions and the crossed beam scattering studies on rare gases with H2O and H2S.

Trifurcated (four-center) hydrogen bond in solid state crystal structure of 5'-amino-5'-deoxyadenosine p-toluenesulfonate

The crystal structure of 5'-amino-5'-deoxyadenosine (5'-Am.dA) p-toluenesulfonate has been determined by X-ray crystallographic methods. It belongs to the orthorhombic space group P2(1)2(1)2(1) with a = 7.754(3)Angstrom, b = 8.065(1)Angstrom and c = 32.481(2)Angstrom. This nucleoside side shows a syn conformation about the glycosyl bond and C2'-endo-C3'-exo puckering for the ribose sugar. The orientation of N5' atom is gauche-trans about the exocyclic C4'-C5' bond. The amino nitrogen N5' forms a trifurcated hydrogen bond with N3, O9T and O4' atoms. Adenine bases form A.A.A triplets through hydrogen bonding between N6, N7 and N1 atoms of symmetry related nucleoside molecules.

Cytotoxicity of half sandwich ruthenium(II) complexes with strong hydrogen bond acceptor ligands and their mechanism of action

Neutral and cationic organometallic ruthenium(II) piano stool complexes of the type [(eta(6)-cymene)R-uCl(X)(Y)] (complexes R1-R8) has been synthesized and characterized. In cationic complexes, X, Y is either a eta(2) phosphorus ligand such as 1,1-bis(diphenylphosphino)methane (DPPM) and 1,2-bis(diphenylphosphino)ethane (DPPE) or partially oxidized ligands such as 1,2-bis(diphenylphosphino)methane monooxide (DPPMO) and 1,2-bis(diphenylphosphino)ethane monooxide (DPPEO) which are strong hydrogen bond acceptors. In neutral complexes. X is chloride and Y is a monodentate phosphorous donor. Complexes with DPPM and DPPMO ligands ([(eta(6)-cymene)Ru(eta(2)-DPPM)Cl]PF6 (R2), [(eta(6)-cymene)Ru(eta(2)-DPPMO)Cl]PF6 (R3), [(eta(6)-cymene)Ru(eta(1)-DPPM)Cl-2] (R5) and [(eta(6)-cymene)Ru(eta(1)-DPPMO)Cl-2] (R6) show good cytotoxicity. Growth inhibition study of several human cancer cell lines by these complexes has been carried out. Mechanistic studies for R5 and R6 show that inhibition of cancer cell growth involves both cell cycle arrest and apoptosis induction. Using an apoptosis PCR array, we identified the sets of antiapoptotic genes that were down regulated and pro-apoptotic genes that were up regulated. These complexes were also found to be potent metastasis inhibitors as they prevented cell invasion through matrigel. The complexes were shown to bind DNA in a non intercalative fashion and cause unwinding of plasmid DNA in cell-free medium by competitive ethidium bromide binding, viscosity measurements, thermal denaturation and gel mobility shift assays.

Peptide mimics for structural features in proteins. Crystal structures of three heptapeptide helices with a C-terminal 6-->1 hydrogen bond.

The crystal structure determination of three heptapeptides containing alpha-aminoisobutyryl (Aib) residues as a means of helix stabilization provides a high-resolution characterization of 6-->1 hydrogen-bonded conformations, reminiscent of helix-terminating structural features in proteins. The crystal parameters for the three peptides, Boc-Val-Aib-X-Aib-Ala-Aib-Y-OMe, where X and Y are Phe, Leu (I), Leu, Phe (II) and Leu, Leu (III) are: (I) space group P1, Z = 1, a = 9.903 A, b = 10.709 A, c = 11.969 A, alpha = 102.94 degrees, beta = 103.41 degrees, gamma = 92.72 degrees, R = 4.55%; (II) space group P21, Z = 2, a = 10.052 A, b = 17.653 A, c = 13.510 A, beta = 108.45 degrees, R = 4.49%; (III) space group P1, Z = 2 (two independent molecules IIIa and IIIb in the asymmetric unit), a = 10.833 A, b = 13.850 A, c = 16.928 A, alpha = 99.77 degrees, beta = 105.90 degrees, gamma = 90.64 degrees, R = 8.54%. In all cases the helices form 3(10)/alpha-helical (or 3(10)helical) structures, with helical columns formed by head-to-tail hydrogen bonding. The helices assemble in an all-parallel motif in crystals I and III and in an antiparallel motif in II. In the four crystallographically characterized molecules, I, II, IIIa and IIIb, Aib(6) adopts a left-handed helical (hL) conformation with positive phi, psi values, resulting in 6-->1 hydrogen-bond formation between Aib(2) CO and Leu(7)/Phe(7) NH groups. In addition a 4-->1 hydrogen bond is seen between Aib(3) CO and Aib(6) NH groups. This pattern of hydrogen bonding is often observed at the C-terminus of helices proteins, with the terminal pi-type turn being formed by four residues adopting the hRhRhRhL conformation.

Influence of the hydrogen bond on the N-H stretching frequencies in amino-acids

The spectra of glycine, its addition compounds and other amino-acids exhibit Raman lines in the region from 3250 cm.−1 to 2500 cm.−1 It has been shown that these lines cannot be assigned to N-H...O stretching vibrations, where the N atom has the covalency of three, but to N+-H...O stretching vibration where the N atom has the covalency of four. Using the data obtained with triglycine sulphate which has the largest number of N+-H...O bonds and whose H bond lengths are known, the correlation curve giving the relation between the N+-H...O stretching frequencies and the corresponding H bond lengths has been drawn. Using this correlation curve, the N+-H...O stretching frequencies appearing inα-glycine,γ-glycine, diglycine hydrochloride, diglycine hydrobromide,l-asparagine monohydrate anddl-alanine have been satisfactorily accounted for on the basis of the known hydrogen bond lengths in these substances.

Raman spectra of single crystals of adipic and sebacic acids and a study of the hydrogen bond vibrations in carboxylic acid dimers

Raman spectra of single crystals of adipic and sebacic acids have been photographed for the first time using λ 2537 excitation. The spectra have been divided into four regions: (a) internal frequencies; (b) summations and overtones; (c) external vibrations; and (d) low-frequency hydrogen bond oscillations. Tentative correlations have been given for all the internal frequencies and summations and overtones. A series of diffuse weak bands observed in the spectra of both these acids in the not, vert, similar2400–2800 cm−1 have been explained as a superposition of O---H frequencies lowered due to hydrogen bond formation over the summations and overtones of fundamentals mainly in the not, vert, similar1000–1500 cm−1 region. Rotatory type of external oscillations of the two formula units of these molecules in their unit cells have been identified at 76, 99, 118 and 165 cm−1 in adipic acid and 66, 95, 117 and 177 cm−1 in the spectrum of sebacic acid. A brief discussion of the low frequency hydrogen bond vibrations in these acids has been made. Making use of the Lippincott—Schroeder potential and assuming a highly anharmonic potential curve for the hydrogen bond, the vibrational frequencies of the bond have been theoretically evaluated. There is very good agreement between these and the experimental values. The results for adipic acid in cm−1 are: 304 (0 → 1), 270 (1 → 2), 241 (2 → 3), 222 (3 → 4) 201 (4 → 5), 183 (5 → 6). In the case of sebacic acid some of the intermediate and higher transitions are absent in the spectrum recorded by the author. From the above data for adipic acid the dissociation energy of the hydrogen bond was evaluated as 5·9 kcal/mole in fair agreement with the values derived from conventional methods.

Evaluation of solute-solvent interactions from solvent blue-shifts of n → π* transitions of C=O, C=S, NO2 and N=N groups: hydrogen bond energies of various donor-acceptor systems

Effects of non-polar, polar and proton-donating solvents on the n → π* transitions of C=O, C=S, NO2 and N=N groups have been investigated. The shifts of the absorption maxima in non-polar and polar solvents have been related to the electrostatic interactions between solute and solvent molecules, by employing the theory of McRAE. In solvents which can donate protons the solvent shifts are mainly determined by solute-solvent hydrogen bonding. Isobestic points have been found in the n → π* bonds of ethylenetrithio-carbonate in heptane-alcohol and heptane-chloroform solvent systems, indicating the existence of equilibria between the hydrogen bonded and the free species of the solute. Among the different proton-donating solvents studied water produces the largest blue-shifts. The blue-shifts in alcohols decrease in the order 2,2,2-trifluoroethanol, methanol, ethanol, isopropanol and t-butanol, the blue-shift in trifluoroethanol being nearly equal to that in water. This trend is exactly opposite to that for the self-association of alcohols. It is suggested that electron-withdrawing groups not merely decrease the extent of self-association of alcohols, but also increase the ability to donate hydrogen bonds. The approximate hydrogen-bond energies for several donor-acceptor systems have been estimated. In a series of aliphatio ketones and nitro compounds studied, the blue-shifts and consequently the hydrogen bond energies decrease with the decrease in the electron-withdrawing power of the alkyl groups. It is felt that electron-withdrawing groups render the chromophores better proton acceptors, and the alcohols better donors. A linear relationship between n → π* transition frequency and the infrared frequency of ethylenetrithiocarbonate has been found. It is concluded that stabilization of the electronic ground states of solute molecules by electrostatic and/or hydrogen-bond interactions determines the solvent shifts.

Potential functions for hydrogen bond interactions - III. Empirical Potential Function for the Peptide N-H...O=C Hydrogen Bond

A careful comparison of the distribution in the (R, θ)-plane of all NH ... O hydrogen bonds with that for bonds between neutral NH and neutral C=O groups indicated that the latter has a larger mean R and a wider range of θ and that the distribution was also broader than for the average case. Therefore, the potential function developed earlier for an average NH ... O hydrogen bond was modified to suit the peptide case. A three-parameter expression of the form {Mathematical expression}, with △ = R - Rmin, was found to be satisfactory. By comparing the theoretically expected distribution in R and θ with observed data (although limited), the best values were found to be p1 = 25, p3 = - 2 and q1 = 1 × 10-3, with Rmin = 2·95 Å and Vmin = - 4·5 kcal/mole. The procedure for obtaining a smooth transition from Vhb to the non-bonded potential Vnb for large R and θ is described, along with a flow chart useful for programming the formulae. Calculated values of ΔH, the enthalpy of formation of the hydrogen bond, using this function are in reasonable agreement with observation. When the atoms involved in the hydrogen bond occur in a five-membered ring as in the sequence[Figure not available: see fulltext.] a different formula for the potential function is needed, which is of the form Vhb = Vmin +p1△2 +q1x2 where x = θ - 50° for θ ≥ 50°, with p1 = 15, q1 = 0·002, Rmin = 2· Å and Vmin = - 2·5 kcal/mole. © 1971 Indian Academy of Sciences.

Potential functions for hydrogen bond interactions - IV. Minimum Energy Conformation of the α-Helical Structure of PoIy-L-Alanine

Making use of the empirical potential functions for peptide NH .. O bonds, developed in this laboratory, the relative stabilities of the rightand left-handed α-helical structures of poly-L-alanine have been investigated, by calculating their conformational energies (V). The value of Vmin of the right-handed helix (αP) is about - 10.4 kcal/mole, and that of the left-handed helix (αM) is about - 9.6 kcal/mole, showing that the former is lower in energy by 0.8 kcal/mole. The helical parameters of the stable conformation of αP are n ∼ 3.6 and h ∼ 1.5 Å. The hydrogen bond of length 2.85 Å and nonlinearity of about 10° adds about 4.0 kcal/ mole to the stabilising energy of the helix in the minimum enregy region. The energy minimum is not sharply defined, but occurs over a long valley, suggesting that a distribution of conformations (φ{symbol}, ψ) of nearly the same energy may occur for the individual residues in a helix. The experimental data of a-helical fibres of poly-L-alanine are in good agreement with the theoretical results for αP. In the case of proteins, the mean values of (φ{symbol}, ψ) for different helices are distributed, but they invariably occur within the contour for V = Vmin + 2 kcal/mole for αP.